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We study an antiferromagnetic spin- 1/2 model with up to third nearest-neighbor couplings on 
the Kagome lattice in the easy-axis limit, and show that its low-energy dynamics are governed by 
a four site XY ring exchange Hamiltonian. Simple "vortex pairing" arguments suggest that the 
model sustains a novel fractionalized phase, which we confirm by exactly solving a modification of 
the Hamiltonian including a further four-site interaction. In this limit, the system is a featureless 
"spin liquid", with gaps to all excitations, in particular: deconfined S' = 1/2 bosonic "spinous" and 
Ising vortices or "visons". We use an Ising duality transformation to express vison correlators as 
non-local strings in terms of the spin operators, and calculate the string correlators using the ground 
state wavefunction of the modified Hamiltonian. Remarkably, this wavefunction is exactly given by 
a kind of Gutzwiller projection of an XY ferromagnet. Finally, we show that the deconfined spin 
liquid state persists over a finite range as the additional four-spin interaction is reduced, and study 
the effect of this reduction on the dynamics of spinous and visons. 



I. INTRODUCTION 

It has now been almost 15 years since P.W. Ander- 
son suggested that two-dimensional (2d) spin one-half 
antiferromagnets might condense into a featureless "spin 
liquid" quantum ground statcEJ. In close analogy with 
the onc-dimcnsional Hciscnbcrg antiferromagnetic chain, 
the 2d spin liquid was posited to support deconfined 
spinon excitations - "particles" carrying s=l/2 in stark 
contrast with the s=l triplet excitations of more famil- 
iar non-magnetic phases such as the spin-Peierls state 
and with the S*^ = 1 magnon excitations of the 2d Neel 
stateH. Early attempts to demonstrate the existence of 
the 2d spin liquid focussed on quantum dimer modelsEI 
motivated directly by Resonating Valence, Bond (RVB) 
ideasQ, slave-Fermion mean field theoriesB and large N 
generalizations^ of the spin models. While the topolog- 
ical character of the spin liquid was mentioned in some 
of these pioneering studiesQ, generally the focus was on 
characterizing the spin liquid by an absence of spin or- 
dering and spatial symmetry breaking. In the past few 
years, it has been emphasized thai the precise way to 
characterize a 2d spin liquid phasdj - as with other 2d 
fractionalized phases - is in terms of "topological order," 
a notion introduced by X.-G. Wen in the context of the 
fractional quantum Hall effectEi. Central to the notion 
of topological order in 2d is the presence of vortexJike 
excitations with long-ranged statistical interactionsElE^I. 
In the simplest 2d spin liquid these point-like excitations 
haveJpeen dubbed "visons" since they carry an Ising or Z2 
fluxllO. Upon transporting a spinon around a vison, the 
spinon's wavefunction acquires a minus sign. A theoreti- 
cal description of this long-ranged statistical interaction 
is most readily incorporated in the context of a gauge 
theory with a discrete Ising symmetry, in which the «i- 
sons carry the Z2 flux and the spinous the Z2 chargeliSl. 
The Z2 gauge theory can be dualized into a vortex repre- 
sentation, wherein the topological order follows from the 



notion of "vortex pairing "El. 

Efforts to identify microscopic spin Hamiltonians that 
might actually exhibit such topologically ordered phases 
have focussed on strongly frustrated 2d s=l/2 antiferro- 
magnets. Due to the "sign problem" these efforts have 
been essentially limited to exact diagonalization stud- 
ies on very small lattices. Nevertheless, such numerics 
do identify a few models which appear to be in a spin 
liquid phase: the Kagome antiferromagnet with near 
neighbor interactionsuHl and a triangular lattice model 
with 4-spin ring exchange termsEj. The importance of 
multi-spin ring exchange processes in driving 2d frac- 
tionalization is palso apparent within the Z2 gauge the- 
ory formulatioiJij. Iri-an important recent development, 
Moessncr and Sondhill3 have compcllingly argued that a 
particular quantum dimer model on the triangular lattice 
is in a featureless liquid phase, closely analogous to a the 
desired "spin liquid" phase of a spin Hamiltonian. 

In this paper we re- visit the s=l/2 Kagome antifer- 
romagnet, in the presence of second and third neighbor 
exhange interactions. By passing to an easy-axis limit of 
this model, substantial analytic and numerical progress 
is possible both in establishing the presence of a frac- 
tionalized spin liquid and of directly analyzing it's topo- 
logical properties. Specifically, in the easy-axis limit we 
map the model exactly onto an XY Hamiltonian con- 
sisting solely of a local 4-spin ring exhange interaction. 
Since the sign of the ring exchange term is "bosonic" 
- opposite to the sign obtained upon cyclically permiit- 
ing four underlying s=l/2 fermions (eg. electrons - 
the Hamiltonian does not suffer from a sign problem and 
so should be amenable to quantum Monte Carlo. Fur- 
thermore, if the two levels of the spin-1/2 on each site 
of the Kagome lattice is reinterpreted as the presence 
or absence of a (quantum) dimer living on a bond of a 
triangular lattice, the model can be reinterpreted as a 
quantum dimer model which is very similar to that con- 
sidered by Moessner and Sondhitil, the distinction being 
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that three, rather than one, dimers emerge from each site. 
This reahzation allows us, to exploit the important work 
of Rokhsar and Kivelsoro who identified an exactly sol- 
uble point of a generalized square lattice quantum dimer 
model. With a similar generalization, our model also 
possesses an exact zero energy wavefunction: an equal 
weight superposition of all allowed spin configurations in 
the low energy singlet sector. We show that this wave- 
function can be viewed as an exact version of the popular 
variational state consisting of |the Gutzwiller projection 
of a superfluid/superconductoiil3. Finally, we are able to 
implement an exact duality transformation which enables 
us to identify the operators which create both the spinon 
excitation and the topological vison excitation. Employ- 
ing the exact wave function, we compute numerically the 
vison 2-point correlation function, and show that it is 
exponentiaUy decaying - the hallmark of a 2d fractional- 
ized phasal3. We thereby demonstrate that the (gapped) 
spinous are genuine deconfined particle-like excitations. 

The paper is organized as follows. In Section II we 
introduce a generalized s = 1/2 Kagome antiferromag- 
net and show how it can be mapped onto a bosonic ring 
model in the easy-axis limit. With a slight further gener- 
alization, we identify an exactly soluble point in Section 
III and obtain an exact spin-liquid ground state wave- 
function. In Section IV we exploit an exact duality trans- 
formation which maps the Kagome spin model onto a Zi 
gauge theory living on the dual lattice to identify the 
spinon and vison excitations. The vison two-point cor- 
relation function is then evaluated numerically using the 
exact wavefunction in Section V, and we demonstrate 
that it is short-ranged thereby directly establishing the 
presence of fractionalization in the spin-liquid ground 
state. Finally, Section VI is devoted to a brief discus- 
sion of the implications of this finding. 



II. MODEL 

We consider a spin one-half Heisenberg antiferromag- 
net on a Kagome lattice with Hamiltonian 



(1) 



Since the Kagome lattice consists of corner sharing trian- 
gles, the nearest neighbor exchange interaction, denoted 
Ji, is strongly frustrating. Here we extend this stan- 
dard near-neighbor model to incudc further neighbor in- 
teractions, JiiJzi which act between pairs of sites on the 
hexagons in the Kagome lattice (Fig. |^). Specifically, two 
spins on the same hexagon separated by 120 degrees are 
coupled via J2, and J3 is the coupling between two spins 
diametrically across from one another on the hexagon. 

Instead of the usual nearest-neighbor Kagome Antifer- 
romagnet (with J2 = "^3 = 0), we specialize instead to the 
case with equal exchange interactions, J\ = J2 = = J ■ 
This generalized Kagome Antiferromagnet can be cast 




FIG. 1: Kagome lattice and interactions. Two primitive vec- 
tors ai , 0,2 are shown, as are the labels 1 ... 4 for the four 
sites on a bow-tie. The ring term is generated both from the 
red (short dashes) and blue (long dashes) virtual exchange 
processes. 



into a simple form: 



TL = J ^ So • So, 

o 



(2) 



where the summation is over all hexagons on the Kagome 
lattice and So = X]i=i Si is the sum of the six spins on 
each hexagop. A similar form was obtained by Palmer 
and Chalkerll3 for a Heisenberg model on the "checker- 
board" lattice, with the Hamiltonian expressed as a sum 
over the total spin living on elementary square plaque- 
ttcs, squared. 

As for the nearest-neighbor model, the generalized 
Kagome antiferromagnet described by Eq. ^ has a non- 
trivial classical limit. There is a thermodynamically large 
set of classical ground states, which includes any config- 
uration for which the classical vector So = for each 
hexagon. The breaking of this degeneracy by quantum 
fluctuations could give rise to "order-by-disorder". For 
the spin- 1/2 case of interest, however, Eq. || is essen- 
tially intractable analytically. To make progress, we re- 
tain SU(2) spins with S — 1/2 on each site, but generalize 
the Hamiltonian to allow for anisotropic exchange inter- 
actions. Specifically, we consider an "easy axis" limit, 
with the exchange interaction along the z— axis in spin 
space larger than in the x — y plane: Jz > J_l. In the ex- 
treme easy-axis limit, one can first analyze the Jz terms 
alone, and then treat the remaining terms as a perturba- 
tion: Ti, = TCo + Hi with. 



and 



T-Cq - Jz'^ (Sq)^ , 



(3) 



(4) 
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where the subtraction of 3 was included for convenience. 
In an eigenbasis oi = S = ±1/2, the Hamiltonian Hq 
describes a classical spin system. The classical ground 
state consists of all spin configurations which have zero 
(z-axis) magnetization on each and every hexagonal pla- 
quette: Sq = 0. There are many such configurations 
(note that unlike the nearest-neighbor model, the gen- 
eralized Kagome antiferromagnet is unfrustrated in the 
easy-axis limit), with a ground state degeneracy that 
grows exponentially with system size, much like other 
fully frustrated classical spin models such as the triangu- 
lar lattice Ising antiferromagnet. The full Hamiltonian, 
H, lifts this huge degeneracy, splitting the classically de- 
generate ground states into a low-energy manifold, still 
characterized, however, by the good quantum numbers 

Some properties of this easy-axis limit are immediately 
evident. For instance, all states in the low-energy mani- 
fold have Sq = for every hexagon, and there is a large 
gap of approximately Jz to states with any non-zero Sq. 
Hence the ground state has in this sense a "spin gap". 
Thus the easy-axis generalized Kagome antiferromagnet 
has no XY spin order, but translational symmetry break- 
ing is not precluded. More subtle aspects of this model 
are less evident. In particular, we would like to ascertain 
the presence or absence of more subtle "topological" or- 
der, and the types of "singlet" (more precisely ~ 0) 
and spinful {S^ ^ 0) excitations. 

To proceed, we treat Tii as a perturbation with J± <C 
Jz, and project back into this low energy manifold of 
degenerate classical ground states with Sq — 0. (This 
procedure is very much analogous to the derivation of the 
Heisenberg model starting from the Hubbard model with 
t ^ U. Indeed, in the language of "hard core bosons" 
in which the boson number corresponds to Sf + 1/2, 
the perturbing Hamiltonian Ti.i describes boson hopping 
amongst a pair of sites on the same hexagon.) Within 
second order degenerate perturbation theory (in Jj_) for 
the low-energy manifold, there are two types of (vir- 
tual) processes which contribute, preserving the vanish- 
ing magnetization on every hexagon. In the first, two 
antiparallel spins within a single hexagon exchange and 
then exchange back again. This "diagonal" process leads 
(within the low-energy manifold) to a simple constant 
energy shift Eo ^ -(9/2)NoJ'i/Jz, where Nq is the to- 
tal number of hexagons. Because this trivial shift does 
not split the extensive degeneracy, we neglect it in what 
follows. More interesting are off-diagonal processes, in 
which two pairs of antiparallel spins on opposite sites of 
a 5-site "bow-tie" plaquette exchange (see Fig. |l|). This 
process involves spins on only four sites, and is an analog 
of electron exchange "ring" moves. One can readily ver- 
ify that such "ring" moves on the bow-tie leave invariant 
the (z-axis) magnetization on every hexagon. 

Up to second order in J±/Jz, within the low-energy 
manifold, the full Kagome Heisenberg antiferromagnet is 



thereby reduced to the form: Tio + Hring with 

'^ring — ~'Aing ^~^('5']^'S'2 S^ S^ +h.C.), (5) 
IX 

where the labels 1 ... 4 denote the four spins at the ends 
of each bowtie as labeled in Fig. |l|. Here the ring ex- 
change interaction Jring = Jj_/Jz, and by assumption 
one has J± ^ Jz- It is noteworthy that in this extreme 
easy-axis limit the frustrated Kagome magnet does not 
have a sign problem, and as such could be profitably at- 
tacked via quantum Monte Carlo. 

III. SOLUBLE SPIN LIQUID 

We now use Tiring to address the nature of the spin- 
gapped state of the easy-axis generalized Kagome anti- 
ferromagnet. Several arguments point to a spin- liquid 
phase which supports fractionalized "spinon" excitations 
which carry spin S"^ = 1/2. Such a fractionalized state 
must also support vortex-like excitations, dubbed "vi- 
sons" , which carry no spin but have a long-ranged statis- 
tical interaction with spinous. 

A first suggestion to this effect comes from viewing 
Tiring as a lattice boson model, and a spin liquid state 
thereby bosonic Mott insulator. Generally, such 

bosonic insulating states can be regarded as quantum- 
mechanical condensates of vorticeslHI. To examine the 
vortex excitations, it is convenient to think of Sf^ as lat- 
tice boson raising and lowering operators. Formally, one 
may then express Sf = e^"^' - fluctuations in the U{1) 
phases (pi (conjugate to Sf) are induced by the constraint 
Sf = ±1/2. It is then illuminating to re-express the 
bosonic ring term as 

Tiring = -2 Jring ^ COs((/)i - (1)2 + (f>3 ~ (6) 

1X1 

Consider now a vortex centered on some site (the "core" ) . 
Classically, for the four sites on the bow-tie surrounding 
the vortex core, = (j/4)27r7V„, where Ny denotes the 
number of vortices (vorticity) on this plaquette. The 
(core) energy of this vortex configuration is proportional 
to 

£;vort = 2Jri„g(l - cos(7V„7r)). (7) 

Notice that plaquettes with an odd number of vortices, 
Ny , cost an energy 4 Jring relative to the even-iVi, plaque- 
ttes. In particular, a single strength vortex is costly, but 
double-strength vortices are cheap. The same conclusion 
can be shown more formally using an exact duality trans- 
formation. 

Typically single strength vortices condense, but one 
can also imagine insulating states which result froni^ 
condensation of composites made from Ny vorticesllil. 
Such insulators are necessarily fractionalized since they 
support deconfined (but gapped) charge excitations with 
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"boson charge" Q = = 1 /Ny . Based on the energetics 
of the ring term which tends to expel single vortices with 
double vortices being energetically cheaper, one expects 
that the insulating state for the Kagome ring model will 
have spin S*^ = 1/2 excitations - if it is fractionalized at 
all. If fractionalized, the "vison" can be understood as 
an unpaired vortex state in the vortex-pair condensate, 
a "dual" analog of a BCS quasiparticle. 

Further evidence that the ground state of this model 
might be fractionalized comes from its formal equiva- 
lence to a particular quantum dimer model. Mapping 
to a dimer model is straightforward since the sites of the 
Kagome lattice can be viewed as the centers of the links 
of a triangular lattice. The two = +(— )l/2 states on 
a site correspond to the presence (or absence) of a dimer 
on the associated link on the triangular lattice. The 
ring term above corresponds directly to the elementary 
quantum-dimer move on the triangiilar lattice considered 
recently by Sondhi and Moessneio. The only difference 
with the standard dimer model is that in this instance 
there are three dimers coming out of every site of the 
triangular lattice instead of the usual one. Sondhi and 
Moessner considered an additional "diagonal" term (see 
below) in the triangular lattice quantum dimer model, 
and argued that the model was in a spin liquid state in 
portions of the phase diagram. Central to their argument 
was an exactly soluble point of tlie model, first exploited 
by Rokhsar and Kivelson (RK)eI in the square lattice 
quantum dimer model. The additional term is diago- 
nal in Sf, and may be written Hni = "4 X^rGcx -^flip('')' 
where 

^flipW= E ri il + <y{-iys]). (8) 

<T=±1 jer=l 

The operator -Pflip(r) is a projection operator onto the 
two flippable states of the bow-tie ring r. This term in 
the Hamiltonian can be combined with Tiring and written 
in the suggestive form: 

4 

Wring + -Hnf = E ^fliP ('^) i - -^""g n + "4 } . (9) 

r j=l 

When U4 — Jiing one can write down exact ground 
state(s) which have the product of 2Sx equalling one on 
all bow-tie rings. One such state is the XY ferromagnet 
with S'J = 1/2 on every site. In the hard-core boson de- 
scription, this corresponds to a superfluid state (albeit an 
unusual one with no zero-point fluctuations). One must 
project back into the subspace in which there arc three 
bosons on every hexagon {Sq = 0), since otherwise this 
state will not be an cigenstate of Ho- (Actually, several 
distinct projections are generally possible, onto different 
sectors disconnected from one another under the action 
of Wring- These give degenerate ground states.) This pro- 
jection of a superfluid wavefunction to obtain a bosonic 
insulating state is analogous to the Gutzwiller projec- 
tions of superconducting wavefunctions, to obtain varia- 
tional states for quantum spin modelsEj, but there is an 



important difference. In the present instance, the con- 
straints (of 3 bosons on every hexagon) commute with 
the Hamiltonian Tiring which hops the bosons, in con- 
trast to the no-double occupancy constraint which does 
not commute with the electron kinetic energy term in 
Hubbard type models. Thus, in our case the wave func- 
tion after projection is still an exact cigenstate of the full 
Hamiltonian. 



IV. DUALITY, VISONS, AND SPINONS 

Before studying this wave function, it is convenient to 
expose the vison degrees of freedom via a duality trans- 
formation. Specifically, we will employ the standard 2-t-l 
dimensional Ising duality which connects a global spin 
model to a Z2 gauge theory ^ith gauge fields living on 
the links of the dual latticella. In our case the global 
spin model is the Kagome model Tiring + Tinf in Eq. |^, 
so that the dual lattice is the "dice" lattice, which can 
conveniently be constructed in terms of two interpene- 
trating honeycomb lattices as depicted in Fig. |[ On the 
operator level, the duality transformation is implemented 
by re-expressing and directly in terms of the dual 
gauge fields, ct^- - a set of Pauli matrices living on the 
links of the dice lattice: 

jieo 

and 

s-^-^E^'^i- (11) 

jl=i 

Here, the first product is taken around an elementary 
four-sided plaquette on the dice lattice which surrounds 
the spin Sf . The second product involves an infinte string 
which connects sites of the Kagome lattice, eminating 
from the site and running off to spatial infinity. For 
every bond of the dual dice lattice which is bisected by 
this string, a factor of afj is present in the product. To 
assure that this definition is independent of the precise 
path taken by the string, requires imposing the constraint 
that the product of af, on all bonds connected to each 
site on the dice lattice is set equal to unity: 

^^ = 114 = 1' (12) 

where here j labels the near- neighbor sites to i. These lo- 
cal Z2 gauge constraints must be imposed on the Hilbert 
space of the dual theory. In the resulting dual gauge 
theory, these constraints are analogous to Coulomb's law 
{\7 ■ E = 0) in conventional electromagnetism. The ne- 
cessity of including the constraints can be simply seen 
by counting degrees of freedom: there are twice as many 
(six) bonds per unit cell on the dice lattice as sites (three) 
per unit cell on the Kagome lattice, hence to maintain the 
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physical Hilbcrt space of the original spins (site variables) 
requires restricting the gauge fields (bond variables). 




FIG. 2: Dice lattice shown as two interpenetrating honey- 
combs, indicated by blue (solid) and red (dashed) lines. A 
blue vison is created geometrically by multiplying 23^ over 
the underlying kagome sites (centers of parallelograms, shown 
by solid dots) through which the "string" indicated passes. In 
the dual variables, this product is given by the product of blue 
gauge fields crfj cut by the string shown. The "blueness" of 
the vison shown owes to the fact that only a single spin is 
contained within the originating blue hexagon. 



The dual Hamiltonian takes the form: 



'Hduai = Alip(?^){-Jri 
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«4}, (13) 



where the products are taken around the two hexagonal 
plaquettes of the dice lattice which surround the given 
Kagome site r. These products measure "magnetic flux" 
(in the dual gauge fields) through hexagons belonging to 
the two honeycomb sublattices. The flip term becomes. 



Aiip(0 



(14) 



where the product is taken over pairs of bonds on the el- 
ementary dice plaquette which both connect to the same 
site, j. 

One can readily verify that the operators which im- 
plement a local gauge transformation, Qi in Eqn. [l^ , 
commute with this dual Hamiltonian. Equivalently, since 



-crfj, the dual Hamiltonian is invariant under 



the general Z2 gauge transformation, 



(15) 



with arbitrary = ±1. Remarkably, though, it turns 
out that this gauge theory actually has an additional set 
of local Z2 symmetries. In particular, it is possible to 



transform the gauge fields living on the blue (or red) 
links separately, and still leave the Hamiltonian invari- 
ant. Eqivalently, one can define local red or blue gauge 
operators which commute with the Hamiltonian: 



(16) 



with the product over red links which eminate from site 
i, and similarly for the blue links. On the dice lattice, 
for each 6-fold coordinate site there corresponds both a 
blue and red local gauge operator, whereas the 3-fold 
coordinated sites are either red or blue. 

The presence of this additional local symmetry can be 
directly traced to the conservation of the magnetization 
Sq on each hexagon of the original Kagome lattice ring 
spin-model (note that this is the conservation of dimer 
number emerging from each site on the equivalent trian- 
gular lattice dimer model). Indeed, upon using Eqn. |l^, 
one can show that for the 6-fold coordinate sites of the 
dice lattice 



exp [iirSf 



OA] 



(17) 



where the center of the hexagon is at site i on the dice 
lattice. The right hand side of this expression can be in- 
terpreted as a Z2 "charge" living on the 6-fold coordinate 
sites of the dice lattice, since it equals the (lattice) diver- 
gence of the Z2 "electric fields" . For the "singlet" sector 
of the theory with Sq = for all hexagons, the right 
hand side is simply minus one. But more generally, this 
expression indicates that hexagons with a non-zero (odd 
integer) value of the globally conserved spin, Sq ~ ±1, 
also carry both a red and a blue Z2 gauge charge. 

This fact allows us to identify both the spinon and vi- 
son excitations in the theory. Specifically, consider start- 
ing in the "singlet" sector of the theory with zero mag- 
netization on every hexagon of the Kagome lattice, and 
flipping a single spin. Since each site of the Kagome lat- 
tice is shared by two hexagons, this creates two hexagons 
each with Sq ~ 1. By adding a small near neighbor 
spin exchange it is possible to hop these two magnetized 
hexagons, and to spatially separate them. As we demon- 
strated above, such magnetized hexagons also carry both 
a red and a blue Z2 charge. Provided the dual gauge the- 
ory is in it's deconfined phase, these magnetized hexagons 
can propagate as independent particles. Since two such 
magnetized hexagons were created when we added spin 
one to the system (by flipping the single spin), each 
magnetized hexagon must carry spin S^ = 1/2, and we 
can thereby identify these excitations as the deconfined 
spinous. 

A 2d spin liquid with deconfined spinous must neces- 
sarily [Siipport topological vortex like excitations - the 
visonsEj. The vison acts as a source of Z2 "flux" for the 
spinon, whose wave function changes sign as it is trans- 
ported around a vison. In the Z2 gauge theory formu- 
lation of 2d fractionalization, the flux of the vison cor- 
responds generally to a plaquette with Hpiaq '^^ ~ ~ -'^ ■ 
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Since the spinons which hop on the 6-fold coordinated 
sites of our dice lattice carry both a red and a blue Zi 
gauge charge, it is clear that this spin liquid phase will 
support two flavors of visons - a red (blue) vison corre- 
sponding to a flux penetrating one red (blue) hexagon of 
the dice lattice. 

Due to the long-ranged statistical interaction between 
visons and spinons, it is not possible to have both exci- 
tations present and freely propagating. In particular, if 
the visons are gapped excitations they will be expelled 
from the ground state and the spinons will be deconfined. 
On the other hand, a proliferation and condensation of 
visons will lead to spinon confinement. Thus, in order to 
establish 2d fractionalization it is adequate to show an 
absence of vison condensation. A useful diagnostic for 
this is the vison 2-point correlation function: 



(18) 



where v denotes a vison creation operator. When this 
correlation function is short-ranged, the visons are not 
condensed, and the system is fractionalized. 

In order to evaluate this correlation function for the 
Kagome spin model, it is necessary to express this vi- 
son 2-point function in terms of the original spin opera- 
tors. To this end, we first note that from the definition 
in Eq. ^ it is apparent that the operator 2S^ creates 
both a red and a blue vison, 2S^ = wfwf, since it intro- 
duces a Z2 flux through the red and blue hexagons of 
the dice lattice which enclose the spin. Since ViVi = 1, 
a single (red) vison (say) can be created by stringing 
together an infintely long product of spin operators S^. 
This "string" starts at the given site of the Kagome lat- 
tice and joins neighboring spins making only ±30*^ turns 
eventually running off to spatial infinity, but otherwise is 
arbitrary. Explicitly, the vison two-point function is then 



k—i 



(19) 



where |0) denotes the ground state, and the product in 
Eq. |l9| is taken, as described above, along some path 
on the Kagome lattice starting at site i and ending at 
site j, containing an even number of sites, and making 
only "±30°" turns left or right. Due to the constraint 
•^hex ~ ^'^^ latter product is path- independent up to 
an overall sign (hence the absolute value in Eq. |l^). We 
also define for convenience the physically interesting spin- 
spin correlator. 



C,, = (0|5fS-|0), 



(20) 



V. CORRELATORS AT THE SOLUBLE POINT 

With V and C defined appropriately in terms of the 
spins, we are now in a position to evaluate them using 
the exact RK wavefunction. Specifically, we consider ex- 
act ground state wavcfunctions (at the RK point) on the 



torus defined by identifying sites connected by the two 
winding vectors Wi = niSi and W2 = n2a2, where ni, 712 
arc positive even integers, and ai, 02 are primitive vectors 
(see Fig. 1^). The degeneracies etc. of such wavcfunctions 
are nearly identical to that discussed by Moessner and 
Sondhi, so we do not go into detail here. We focus on 
the wavefunction |0), which has been projected onto a 
single topologically connected sector. 

The expectation values of interest can be evaluated 
stochastically using a classical "infintc temperature" 
Monte Carlo algorithm, which "random walks" through 
the various components of the wavefunction. Our numer- 
ical results for a torus with 711 = 712 = 20 are shown in 
Fig. H plotting In V, In C versus distance. Apart from a 
saturation due primarily to roimd-off error, both corre- 
lators clearly display exponential decay. In Cy , In Vij ~ 
— \ri ~ rj\/£^ with apparently the same correlation length 

Short-range exponentially decaying correlations in Cij 
establish the absence of spin order, but do not preclude 
broken translational symmetries such as plaquette or 
bond order. The exponential decay of the vison corre- 
lator V, however, implies that the phase is necessarily 
fractionalized with deconfined Sz = 1/2 spinon excita- 
tions, regardless of the presence or absence of broken 
translational symmetries. The exponential spatial de- 
cay of Vij is suggestive of a vison gap, i.e. the existence 
of a minimum energy required to excite a vison. A vi- 
son gap, however, strictly speaking requires exponential 
decay of the vison correlator in imaginary time, not in 
space at equal time. Conceivably, the latter condition 
could occur in the absence of a vison gap, provided that 
these visons were localized. Given the peculiarities of 
the present model, we desire a direct argument for a vi- 
son gap. Fortunately, such an argument can also be made 
using the properties of the RK point. Our construction 
closely follows and only slightly refines an argument used 
in Ref. |l^. We consider a family of models defined by sin- 
gling out a single "central" triangle of the Kagome lattice, 
and the associated three bow-tie ring moves centered on 
the three sites of this triangle. For these three ring terms, 
we let U4 = ^I'iiig vary independently of U4 — Jring for all 
other bow-ties. Independent of the choice of J^'ing — 0, 
the original, translationally- invariant |0) state remains a 
zero energy ground-state wavefunction. However, for the 
special case J^'j^g = 0, we can find an additional exact 
zero-energy eigenstate, 



(21) 



where Vi is the (string) vison creation operator emerging 
from any of the sites of the central triangle. Thus for 
"^r'ing = 0, there is no vison gap. We interpret this re- 
sult to mean that by reducing the ring couplings on the 
central triangle, a vison has been bound to this triangle, 
with a binding energy that exactly equals its gap in the 
bulk. To test this hypothesis, we calculate the first-order 
energy shift as J^'j^g is increased from zero to positive 
values using perturbation theory. To leading order, one 
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FIG. 3: Correlators for the (20, 20) Kagome torus. The hori- 
zontal axis min|£i — Xj \ is the shortest distance between sites 
i and j on the torus, in units of the inter-site distance. Large 
circles with error bars show the logarithm of Vij , while small 
circles represent ln(Ci-,). Each small circle represents an av- 
erage over all pairs of sites with a fixed separation, while the 
large circles with error bars represent the distribution of string 
correlators of pairs of sites connected by a horizontal string. 
Both correlators clearly decay exponentially, with apparently 
the same exponent, albeit with larger numerical errors for the 
string correlator, due presumably to the fact that this data is 
less spatially averaged. 



finds that 



= {v\n'\v) 

= {0\V^H'V^\0) 

= 6j;i„g(o|Pflip|o), 



(22) 



XY ring model, and moreover shown that, with the addi- 
tion of the M4 interaction, this model is in a topologically- 
ordered phase at the RK point where = Jiing. At this 
point, the model is a true "short-range" spin liquid, in- 
sofar as there is evidently no order or broken symmetry, 
and all excitations are gapped. But more importantly 
by computing the vison 2-point function we explicitly 
demonstrate that this spin-liquid phase is fractionalizcd, 
and supports = 1/2 spinon excitations whose gap is 
0{Jz). Since the visons also are gapped - with a gap of 
O(Jring) = 0{J]_/Jz) - this spin Hquid state is resihcnt. 

In particular, we may consider a variety of perturba- 
tions away from the special soluble model. Followiag 
similar arguments to those of Moessner and SondhiEJ, 
the spin liquid state remains the ground state for A < 
Ui/Jring < 1, where A is some unknown dimensionless 
number. If ^ < 0, then the fractionalized phase persists 
to the pure ring model, but it is possible that an inter- 
mediate phase (or phases) intervene (s) such that A > 0. 
Quantum Monte Carlo simulations could be useful in de- 
ciding this question. 

Perhaps more novel perturbations consist in devia- 
tions from the condition Ji = J2 = J3 imposed initially. 
Again, the presence of a complete gap in the spectrum 
rules out the destruction of the spin liquid by these per- 
turbations (provided they are weak). It is interesting 
to consider the effect of small changes in Ji, in particu- 
lar Jiz = Jz + 5Jiz, Ji_L = + SJi^. Viewing these 
deviations as perturbations, the change in the easy-axis 
coupling can be rewritten as 



(24) 



where H.' is the sum of ring terms for the three sites of the 
triangle. To obtain the final result, we use ViSjVi = ±S'J, 
where the minus (plus) sign obtains if j is (is not) on 
the string. The energy of the no-vison state Eq ^ Q for 
all Jring- Hence the gap in this approximation is pro- 
portional to Jj'jjjg multiplied by the probability that any 
given bowtie is flippablc. The latter probability is deter- 
mined directly by the classical Montc-Carlo procedure, 
and wc find (0|Pflip|0) « 0.257, hence 



E„ 



1.54 Jj.;jjg. 



(23) 



Note that this is only the first-order approximation to 
the vison gap. The naive extension J^i^g Jrins gives 
a reasonable extrapolation (this is used in Ref. |l9|), but 
there is no obvious reason to expect it to be exact. 



VI. DISCUSSION 

We have demonstrated the equivalence of the general- 
ized Kagome antifcrromagnct in the easy-axis limit to an 



where the sum in the first line is taken over nearest- 
neighbor site of the Kagome lattice, and is equivalent in 
the second line to a sum over hexagons that are nearest- 
neighbors within cither the red or blue honeycomb sub- 
lattice of the dual dice lattice. Remarkably, the latter 
form, Eq. ^ corresponds to a vison hopping or kinetic 
energy term. Because the visons are already gapped, this 
clearly will not destabilize the ground state provided the 
kinetic energy gain remains small relative to the vison 
gap, i.e. SJiz < Jring- 

The change in the in-plane exchange can be written 



6Ha 



(25) 



For any given bond on the Kagome lattice, the associated 
term in Eq. ^ raises Sq of one hexagon by +1 and lowers 
of another by —1. Clearly, acting upon the ground 
state, this takes the system outside the low-energy man- 
ifold of Sq — 0. Hence due to the large spin gap, it 
generates only weak second-order virtual processes that 
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renormalize couplings of the effective ring model. How- 
ever, its effects are more interesting on some of the ex- 
cited states. In particular, for a single spinon excitation, 
one has for SJi± = a single magnetized hexagon with 
Sq = ±1. The spinon is completely localized, and there 
is an associated degeneracy of these excited states reflect- 
ing translational invariance due to the arbitrariness of 
which hexagon is magnetized. For Jij_ ^ 0, this degener- 
acy is lifted, since Eg. p5| allows the magnetized hexagon 
to move. Thus Eq. ^ has the effect of giving the spinous 
some kinetic energy, and the associated states broaden 
into a band. Indeed, it is possible to formally rewrite 
Eq. explicitly as a spinon hopping term. As above, 
because of the existing spinon gap, the ground state is 
expected to be stable to this perturbation for SJi± ^ J^. 

We conclude with a comparison of our results to some 
related discussions in the literature. An interesting as- 
pect of the spinous in the generalized Kagomc antifcr- 
romagnet we have considered is that they arc bosonic. 
Despite the close relation of the topologically ordered 
state described here to a Gutzwiller-projected superfluid, 
this is in contrast to what is obtained by such-|projec- 
tions on SU(2)-invariant superconducting stategl3, as are 
naturally suggcS|tpd by work arising from various slave- 
fermion theoricda. As in|-.our work, the large A^ ap- 
proaches to the spin liquidD also find bosonic spinous. 

One of the most intriguing aspects of the numer- 
ical results on pUjie spin-1/2 nearest-neighbor Kagome 
antiferromagneltJ is the proliferation of a large number 
of very low-energy singlet excitations. Our approach does 
not shed too much light on this phenomenon, since the 



spin liquid ground state found here is in fact fully gapped. 
It is, however, true that in the easy-axis limit consid- 
ered above the visons (which are "singlets" under U(l) 
rotations about the 5^ axis) have a much smaller gap 
("/ring = J±/Jz *C Jz) than the spinous. Moreover, it is 
natural to expect in our effective ring model that as the 
ratio 1*4/ Jring is decreased, some confinement transition 
should occur. At such a confinement transition, provided 
it is second order, the vison gap must vanish. Should the 
pure ring model lie near to this critical point, one would 
indeed expect a large number of low-lying singlet excita- 
tions, which on the deconfined side of the critical point 
are understood as weakly gapped visons. 

In light of these results it should be interesting to look 
at other bosonic ring models, well away from the inte- 
grable RK point. For instance, the properties of the pure 
XY ring model, Eq. ^ defined on four-site plaquettes on 
diverse lattices (Kagomc, square, triangular...) arc very 
poorly understood. With the insight that such terms 
strongly favor vortex pairing, these seem excellent can- 
didate models that might exhibit quantum number frac- 
tionalization. A Asiriety of numericalES and novel ana- 
lytical techniquesEll might profitably be applied to these 
systems. 
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